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ABSTRACT 

We present a complete spectroscopic survey of 2414 2MASS-selected blue horizontal branch (BHB) 
candidates selected over 4300 deg^ of the sky. We identify 655 BHB stars in this non-kincmatically 
selected sample. We calculate the luminosity function of field BHB stars and find evidence for very few 
hot BHB stars in the field. The BHB stars located at a distance from the Galactic plane \Z\ < A kpc 
trace what is clearly a metal-weak thick disk population, with a mean metallicity of [Fe/H] = —1.7, a 
rotation velocity gradient of dvrot/d\Z\ = —28 ± 3.4 km s"'^ in the region \Z\ < 6 kpc, and a density 
scale height of hz = 1-26 ± 0.1 kpc. The BHB stars located at 5 < |Z| < 9 kpc are a predominantly 
inner-halo population, with a mean metallicity of [Fe/H] — —2.0 and a mean Galactic rotation of 
—4 ± 31 km s^^. We infer the density of halo and thick disk BHB stars is 104 ± 37 kpc~'^ near the 
Sun, and the relative normalization of halo to thick-disk BHB stars is 4 ± 1% near the Sun. 
Subject headings: stars: early types — stars: horizontal-branch — Galaxy: stellar content — Galaxy: 
halo 



1. INTRODUCTION 

Theoretical simulations show that the remnants 
of hierarchical galaxy formation in the Milky Way 
should sti ll be visible a s star streams in th e stel- 
lar halo (I Johnston et all 119961: iHarding et all l2001l: 



Abadi et al.ll2003l: iBullock fc Johnstonll2005l: iFont et all 
200611 . Star counts and color-magnitude diagrams 



have proven very effective in identifying structures in 
the halo, includ ing the Sa gittarius stream wrapping 
around the sky (jMajewskTet al.. ,2003ll and overdensi- 
ties in Monoceros (iNewberg et al.l 120021 : llbata et ah 
20031 lYannv et all I2003D . Canis Major (iMartii] 



20041) . Tria ngulum- Andromeda (iRocha-Pinto et al 
200l. Virgo (iDuffau et all [20061 : |Vivas & Zinn' '2006'; 
Newberg et al.l 120071) . and e l sewhe re (^ elokur ov et al.. 
20061 iGrillmair k Dionatoi l200l iGrillmaiil 12006^ 



Stellar spectroscopy opens up the possibility of finding 
structures i n velocity ( s uch a s the Sagittarius dwarf 
galaxy, e.g. llbata et al.l Il994j) . in metallicity, and in 
distance. The major difficulty in mapping the stellar 
halo is finding tracer stars that are luminous enough to 
observe at great depths yet common enough to densely 
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samp le the halo. 

In .Brown et al.l (|2003l hereafter Paper I), we intro- 
duced the Century Survey Galactic Halo Project, a pho- 
tometric and spectroscopic survey from which we se- 
lected Blue Horizontal-Branch (BHB) stars as probes 
of the Milky Way halo. BHB stars are evolved, he- 
lium core-burning s tars that serve as e xcellent "stan- 
dard candles." In iBrown et al.l (j2005l . hereafter Pa- 
per II), we explored the Two Micron All Sky Survey 
i|Skrutskie eFali l2006l 2MASS) and th e Sloan Digital 
Skv Survev (lAdelman-McCarthv et all l2006l . SDSS) as 
the basis for a large spectroscopic survey of BHB stars. 
In Paper II we calculated the first field BHB luminosity 
function, and concluded that field BHB stars are consis- 
tent with populations seen in most globular clusters, but 
inconsistent with globular clusters that have substantial 
extended BHBs. 

Here we describe a complete, non-kinematically se- 
lected sample of BHB stars cov ering 10% of the en tire 
sky. Our survey is inspired by iBrown et al.l (|2004D , in 
which we photometrically selected BHB candidates from 
the completed 2MASS catalog. We have now obtained 
spectroscopy for 2414 2MASS-selected BHB candidates, 
allowing us to measure velocities and metallicities for 
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Fig. 1. — AitofT sky map in Galactic coordinates, showing the 
number counts of 2MASS-selected BHB candidates. Pixels are 1 
deg^ in size. Solid thick lines indicate our survey regions. White 
regions in the disk are regions of high reddening that are excluded. 



stars to a depth of 8 kpc over a 4300 deg^ region. Such 
a large-area survey is necessary to unambiguously iden- 
tify halo structure: theoretical simulations predict that 
star streams cover hundreds of square degree s on the sky 
(|Bullock fc Johnston|[2005HFont et al.ll2006f) . 

Pre v ious spectroscopic surv e ys of fi eld BHB stars (iPier 
1983t ISom mer-Larsen et alJ 1981 lArnold &: Gilmora 
1992t iKin man et al. 19 9J; iWilhelm et al.l Il999bl : 
Kinman et al.| 120041: IClewlev et al.l [200l 120051 : 
Kinman et al. 1 I2007fl have identified BHB stars over 
relatively small fractions of the sky (10^ - 10'^ deg^) 
compared to the Century Survey Galactic Halo Project. 
The exception is the sample of 1170 BHB stars observed 
by the SDSS as mis-identified quasars or as fill er 
objects in low density regions (iSirko et al.l [2004al lbl: 
IClewlev fc Kinmanl l2006l ) . In comparison, our spectro- 
scopic survey of BHB stars is cleanly selected and is 
100% complete within the selection limits. 

Although our large-area spectroscopic survey is a rich 
source for general studies of the thick-disk and inner-halo 
populations, here we focus our attention on the proper- 
ties of the BHB stars. Our goal is to lay the groundwork 
for a structure analysis to be presented in later paper (in 
preparation). In §2 we describe spectroscopic observa- 
tions of stars in the new 4300 deg^ region, and discuss 
our selection efficiency for BHB stars. In §3 we discuss 
the global properties of the sample, including the mean 
galactic rotation and metallicity distribution of the stars. 
In §4 we calculate the luminosity function of our clean 
sample of field BHB stars. In §5 we fit for the density 
distribution of the BHB stars and determine the rela- 
tive normalization of thick disk to inner halo stars. We 
conclude in §6. 

2. DATA 

2.1. Target Selection 

Following iBrown et al.l ((2004^ , we select candidate 
BH B stars by color from the 2MASS point source cata- 
log (|Skrutskie et al.l[2006l) ^ with -0.2 < ( J - H)o < 0.1 
and —0.1 < {H — K)o < 0.1. We use de- reddened colors 
and magnitudes (jSchlegel et al.l [19981 ) to create a clean 
sample. The color selection is designed to provide a rel- 
atively high selection efficiency (~40%) for BHB stars, 
but a reduced completeness for BHB stars. Comparison 

^ Available at http: //www. ipac . caltech.edu/2mass/] 
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Fig. 2. — The four BHB classification techniques applied to 
our sample: (a) the modi fied [Kinman et al . (1994) method, (b) 
the IWilhelm etall Ill999al) method, (c) the Clewley eTaT] 11200% 
i5o.i5-Color method, and (d) the[Clewlcv et al. (2002) Scale width- 
Shape method. We consider objects satisfying 3 or more of the 4 
techniques BHB stars [solid symbols). 



with the Paper I sample suggests that th e color selec- 
tion s amples ^67% of the BHB population ([Brown et al.l 
[200l . 

We select BHB candidates in the magnitude range 
12.5 < Jo < 15.5. Our goal is to sample BHB stars 
as distant as possible, yet at J = 15.5 the uncertainty in 
{J — H)o exceeds ±0.1 and thus there is little point in go- 
ing fainter than J = 15.5 with 2MASS. We set J = 12.5 
as our bright limit to avoid thin disk contamination; a 
typical BHB star with My = +0.6 is 2 kpc distant at 
J = 12.5. 

Figure [l] an Aitoff sky map plotted in Galactic co- 
ordinates, shows our survey region. The greyscale indi- 
cates the number counts of 2MASS BHB candidates in 
the magnitude range 12.5 < Jq < 15.5. Our survey re- 
gion includes the north Galactic cap opposite the bulge 
90° < ; < 270°, 60° < 5 < 90° plus an extension to 
b > 35° at 90° < I < 135°. In the south, our survey sam- 
ples a similar region bounded by & < —35°, I > 90°, and 
Dec > —10°. The survey areas cover 2136 deg^ in the 
north Galactic hemisphere and 2190 deg^ in the south 
Galactic hemisphere. 

There are 2414 BHB candidates in our survey region. 
The average surface density of 2MASS-selected BHB can- 
didates is 0.56 deg~^. Our survey is 100% complete and 
can identify stars moving at any radial velocity. 

2.2. Spectroscopic Observations 

Spectroscopic obser vations were obtai ned with the 
FAST spectrograph fFa bricant et al.lll998f ) on the Whip- 
ple 1.5m Tillinghast telescope. Observations were ob- 
tained over the course of 48 nights in 2004 and 2005. 
The spectrograph was operated with a 600 line mni^^ 
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Fig. 3. — Fraction of all A-type stars and BHB stars in our 
sample as a function of apparent J magnitude. 



Fig. 4. — Fraction of all A-type stars and BHB stars in our 
sample as a function of Galactic latitude \b\. 



grating and a 2 arcsec slit, providing spectral resolution 
of 2.3 A and wavelength coverage from 3450 to 5450 A. 
Exposure times were chosen to yield a typical signal-to- 
noise S/N = 30 in the continuum. 

Paper I contains details of the data reduction. We use 
the spectra to measure radial velocities, spectral types, 
metallicities, effective temperatures, and surface gravi- 
ties of the 2414 BHB candidates. During the course of 
this survey we re-observed 30 objects from Paper I. The 
scatter in the spectroscopic measurements of the same 
objects provides us with a direct measurement of the un- 
certainties: ±16 km in velocity, ±1.2 sub-types in 
spectral classification, ±0.4 dex in [Fe/H], ±400 K in 
effective temperature, and ±0.3 dex in surface gravity. 

2.3. BHB Classification 

The major difficulty in using BHB stars as probes 
of Galactic structure is the need to distinguish reli- 
ably between low surface-gravity BHB stars and higher 
surface-gravity A-type dwarfs and blue stragglers. Al- 
though investigators once thought blue stragglers were 
a minor component of the ha lo population, recent 
studies (i Norris fc Hawkins! 119911: iKinman et all 11994 
Preston ct all 119941 : IWilhelm et al.(ll999bl: iBrown et all 
20031 1200|) demonstrate that a surprisingly large frac- 
tion of faint stars in the color range associated with BHB 
stars are indeed high- gravity stars , many of which are 
blue stragglers (,Preston fc SnedenI [200n[ iCarnev et al.l 
[20051) . 

Our classification of BHB stars is identical to the ap- 
proach described in Paper I. In brief, we apply t he tech - 
niques oflKinman e t al. (199j), IWilhelm et all ()1999af ). 
and IClewlev et 31^20021 . i200l to identify low surface- 
gravity BHB stars. Figure [2] displays the results for our 
sample. We consider objects that satisfy 3 or more of the 
4 classification techniques as BHB stars (solid symbols in 
Figure [2]); we identify 779 probable BHB stars. 

We expect halo stars to be largely a metal-poor pop- 
ulation (e.g., Paper II), yet 124 (16%) of the BHB stars 
are relatively metal-rich [Fe/H] > —0.6. Curiously, the 
BHB stars with [Fe/H] > —0.6 are systematically 0.06 
mag bluer in (B — V)a, or 600 K hotter, than the more 
metal-poor BHB stars. Hot BHB stars have weak Ca ii 
K making metallicity measurements difficult. Further- 



more, BHB and main-sequence A stars have similar sur- 
face gravities at ~ 10* K, making classification difficult. 
Thus we consider the [Fe/H] > —0.6 BHB stars suspect 
and mark them as BHB / A stars in Fig. [51 

Because our goal is to create a clean sample of non 
thin-disk BHB stars, we consider the 655 probable BHB 
stars with [Fe/H] < —0.6 as our "clean" s ample of BHB 
stars. This is consistent with lChiba fc Bee rs (2000), who 
use [Fe/H] < —0.6 to select thick-disk and halo stars by 
metallicity. We use the clean sample of BHB stars for 
the following analyses. 

2.4. Sample Selection Efficiency 

Our net selection efficiency for BHB stars is 27% (655 
of 2414), and is a function of both depth and Galactic 
latitude. Figure [3| plots the fraction of all stars of spec- 
tral type A and the fraction of BHB stars as a function 
of apparent J magnitude in our sample. Our color se- 
lection is 90% efficient for selecting A-type stars in the 
interval 12.5 < J < 14.5, but the efficiency plummets 
to 40% at J = 15.5 due to increased photometric errors. 
Interestingly, the relative fraction of A-type stars identi- 
fied as BHB stars increases from 30% at J = 13 to 40% 
at J = 15. The increasing percentage of BHB stars re- 
flects the relative fraction of different stellar populations 
at different depths in the halo. This behavior is best il- 
lustrated in Figure [4l which displays the fraction of all 
A-types and BHB stars in our sample as a function of 
Galactic latitude. BHB stars comprise 20% of the en- 
tire sample with 35° < |6| < 45° and 50% of the entire 
sample with 80° < |6| < 90°. 

The non A-type stars in our sample are mostly early F- 
type stars scattering into our color-selection region, plus 
a small number of hot subdwarfs and white dwarfs with 
blue colors. We classify 46 objects (2% of the sample) as 
subdwarfs and 21 objects (1% of the sample) as DA white 
dwarfs. One object, CHSS 3842 (2MASS J010324.54- 
063210.5), is a hot DB white dwarf. 

We present the photometric and spectroscopic param- 
eters for all 2414 stars in Appendix A. 

2.5. Variable Stars 

We don't expect to find RR Lyrae variables in our sur- 
vey because our color selection targets stars bluewards 
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Fig. 5.— 2MASS J002334.02+065647.6, possibly an accreting 
white dwarf. 



of the horizontal branch instabihty strip. That said, our 
spectroscopy shows that 20% of the sample is composed 
of redder, F-type stars, and our survey is well matched 
to existing variability surveys such as the All Sky Auto- 
mated Survey (ASAS,|Poimanski 200 ^) and the Nort hern 
Sky Variability Survey fNSVS. Wozniak et al.ll200l . 

We match our entire target list to ASAS and find three 
variables: two RR Lyraes (CHSS 3168 and CHSS 3704) 
and one eclipsing binary (CHSS 3341). The NSVS is 
better matched to our survey region and survey depth. 
We find 1519 stars with NSVS photometry, of which six 
have RMS photometric scatter greater than 3 times their 
median photometric error. Of the six possible NSVS vari- 
ables, two are RR Lyraes (CHSS 1983 and CHSS 2983) 
and four show no clear periodicity. Thus our survey con- 
tains a total of four known RR Lyraes, all of which fall 
in the reddest 15% of the sample. The four RR Lyraes 
are not part of the BHB sample. 

2.6. A Possible Accreting White Dwarf 

The most unusual object in the sample is CHSS 3134 
(2MASS J002334.02-H065647.6), possibly an accreting 
white dwarf. The object is ident ified as a DA white dwarf 
in the iBerger fc FringantI 1)19801 ) catalog of blue objects. 
Our higher-resolution spectrum (see Figure [5|), however, 
shows hydrogen Balmer emission lines in the cores of the 
absorption lines from H/3 to HID. Be stars have simi- 
lar spectra, but Be stars usually have Balmer em ission 
in only Ha and H/3 (e.g., iBragg fc Kenvonll2002f ). The 
presence of strong Hei 4471 indicates the star is hot, 
with Toff ?J 25, 000 K, and the broad Balmer absorption 
lines indicate the star has high surface gravity. Based 
on the observed spectrum (Figure [5]), this star is possi- 
bly a white dwarf accreting matter at low rates from a 
close binary companion. A white dwarf can dominate 
the spectrum if mass transfer has reduced its donor star 
to almost nothing. Alternatively, this star could be a 
compact binary in which the white dwarf is illuminating 
a low-mass companion. Further spectroscopic follow-up 
is needed to establish the exact nature of this unusual 
system. 

3. GLOBAL PROPERTIES 

3.1. Disk, Halo Models 

Given the location and depth of our sample, we ex- 
pect our BHB stars to sample both the thick-disk and 
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Fig. 6. — Relative contribution of thin disk, thiclt disk, 
and halo populati ons for two representative star count models 
llSiegel et al.ll200g) calculated for our survey volume. Solid line: 

Zo,thin = 230 pC, Ro^thin = 2 kpc, Pthick = 10%, ^0, thick = 600 

pc, Ro,thiGk = 4 kpc, phalo = 0.15%, c/ahalo = 0.5, r^„^,°. 
Dashed line: Zo^thin = 230 pc, Ro,thin = 2 kpc, pthick = 10%, 
Zo,thick = 600 pc, Ro,thick = 3 kpc, phaio = 1%, c/a^aio = 0.7, 
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halo populati ons. The halo may n ot be a single entity, 
and the recent ICarollo et al.l (j2007[ ) analysis of the abun- 
dances and kinematics of ~20,000 stars from SDSS paints 
a picture of 1) a fiattened inner halo with little or no ro- 
tation and with peak metallicity around [Fe/H] = —1.6, 
and 2) a more spherical outer halo population that is 
strongly counter-rotating and with a peak metallicity of 
around [Fe/H] — —2.0. Our sample of BHB stars reaches 
heliocentric distances up to 10 kpc. Thus, in this pic- 
ture, the majority of our halo stars are associated with 
the "inner-halo" component of the halo. 

To provide context for discussing the properties of our 
sample, we begin by investigating the expected contri- 
bution of the thick-disk and halo populations. Figure [6] 
shows star-count predictions from lSiegel et al.l (|2003) for 
the relative contribution of the thin disk, the thick disk, 
and the halo in our survey volume. The two sets of lines 
(solid and dashed) illustrate the range allowed by their 
best-fit parameters (see their Table 6). 

Figure [6] shows that while the relative contribution of 
the thick-disk and halo populations in our survey volume 
is uncertain, the thin-disk population should be negligi- 
ble. The star-count models suggest that the thick disk 
should dominate our survey for \Z\ < 4 kpc, while the 
inner-halo should dominate for |Z| > 6 kpc. Given the 
uncertainties in the normalizations, however, we simply 
conclude that our survey contains a mix of thick-disk and 
halo stars, and that the halo contribution increases with 
distance from the plane. 

3.2. Radial Velocities 

Figure [7] displays our heliocentric radial velocities, cor- 
recte d for Solar moti on relative to the local standard of 
rest (|Hogg et al.|[20ba) . as a function of spectral type. 
The large group of stars near AO are the BHB stars 
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Fig. 7. — Spectral types and heliocentric radial velocities with 
respect to the local standard of rest. Solid symbols mark the BHB 
stars. Errorbar indicates the average uncertainty of the measure- 
ments. 
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Fig. 8. — Distribution of metallicity [Fe/H] and heliocentric 
velocity corrected to the Local Standard of Rest f0,LSR fof our 
sample of stars. BHB and BHB/A stars are plotted on the left; 
non-BHB stars are plotted on the right. Errorbar indicates the 
average uncertainty of the measurements. 



(solid symbols). The other B-type objects are possibly 
hot BHB stars, blue stragglers, or run-away B stars. The 
F- and late A-type stars exhibit a smaller velocity dis- 
persion than the BHB stars, consistent with their being 
mostly nearby stars located in the disk. All of the high 
velocity stars are probably halo stars. Because our sam- 
ple covers a large area of the sky, we must remove the 
effects of Galactic rotation before calculating the line-of- 
sight velocity dispersion of the halo stars. We calculate 
velocity dispersion and Galactic rotation below. 

3.3. Proper Motions 

It would be very interesting to know the full space 
motions of our stars. We search existing proper mo- 
tion catalogs and find 703 matches wit h the US Naval 
Observatory CCD Astrograph Catalog f Zacharias et al.l 
[2004, UCAC2), 955 match es wit h the US Naval Obser- 
vatory Bl Catalog (jMonet et all l2003l USNOBl), and 
2414 matches with the Guide Star Catalog 2.3 (GSC2.3, 
B. McLean, 2005 private communication). We proceed 
cautiously, however, because our stars are relatively dis- 
tant, at 2 < d < 10 kpc, and the reported proper motions 
are typically quite small, ~10 mas yr~^. 

We compare proper motions between the UCAC2, US- 
NOBl, and GSC2.3 catalogs. The UCAC2 and GSC2.3 
proper motions correlate well, but USNOBl proper mo- 
tions are systematically discrepant from the other two 
catalogs for proper motions less than 10 mas yr~^. The 
UC AC2 appears the most reliable of the three catalogs 
(^see lMink et al.ll2004f) . However, even the UCAC2 may 
contain systematic errors on the scale of degrees (N. 
Zacharias, private communication, 2005), making com- 
parison of stars in different parts of the sky problematic. 
The dispersion of the proper motions between the three 
catalogs is ±7 mas yr~^; we consider this estimate a good 



measure of the accuracy of the catalogs over large areas 
of sky. 

If we restrict ourselves to proper motions with > 3cr 
confidence, the number of stars with proper motions 
greater than 20 mas yr~^ is approximately 10% of the 
catalog matches in all three catalogs. We estimate tan- 
gential velocities for these stars by combining the proper 
motions with our distance estimates (see below). The 
formal uncertainty of the tangential velocities is approx- 
imately ±300 km s"""^. Because the uncertainty greatly 
exceeds expected stellar velocities, these tangential ve- 
locity estimates based on the best proper motions are in 
practice useless for our analyses. Thus we ignore proper 
motion in the remainder of our paper. 

3.4. Metallicities 

The strongest indicator of metallicity in our A-star 
spectra is the 3933 A Ca ii K line. We estimate stellar 
metallicities based on Ca ii as described in Paper I. In 
brief, we use three different techniques: the spectral line 
indices of Beers et al. (1999), the equivalent width of Ca 
II K plus a chi-square comparison betwe en metallic-line 
regions in synthetic and observed spectra (jWilhelm et al.l 
|T999a) , an d a Nelder-Mead algor i thm that fits the entir e 
spectrum (jNelder fc MeadI 119651 : lAllende Prietol 120031) . 
The final metallicities are the average of the three tech- 
niques and have formal uncertainties of ±0.25 dex. As 
mentioned above, we re-observed 30 objects from Paper 
I and find that their metallicity determinations have a 
±0.4 dex RMS scatter. 

Figure [5] plots the metallicities and velocities of the 
BHB and BHB/A stars {left panel) and all the other non- 
BHB stars (right panel) in our sample. The extra stars 
at [Fe/H] = —3 and are artifacts of our method; our 
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measurements are restricted to —3 < [Fe/H] < 0. 

Examination of Figure [8] reveals that the BHB stars 
are more metal-poor than the other stars in our sam- 
ple. Excluding the stars at the [Fe/H] limits, the BHB 
stars have a median [Fe/H] — —1.7. By comparison, the 
non-BHB stars in our sample are more metal-rich, with a 
median [Fe/H] = —0.8. The distributions of metallicity 
and velocity are consistent with the BHB stars consti- 
tuting a largely halo population and the non-BHB stars 
constituting a largely thick-disk population. 

3.5. Distances 

BHB stars are approximate standard candles with lu- 
minosities dependent on effective temperature (color) 
and on metallicity. We estimate {B — V)o, which we la- 
bel BVO, for our BHB stars using 2MASS photometry 
and Balmer line strengths (and SDSS photometry, where 
available), as described in Paper I. We then calculate lu- 
minosities for our BHB stars b y ada pting the My (BHB) 
relation from IClewlev et al.l (|2004D . This relation as- 
sumes the Hipparcos-derived zero point, Mv(RR) = 
0.77±0.13 at [Fe/H] = -1.60 (|Cxould fc PopowskilllQQSh . 
a My-metallicity slope 0.214 ± 0.047 based on RR Lyrae 



stars in the Large Magellanic Cloud (Clcmcnti ni et al 



2003) ■ and a cubic relation in {B — V)q (|Preston et al 



1991) to provide a temperature correction. A detailed 



comparison of the luminosity function of field BHB stars 
and globular cluster BHB stars in Paper II revealed a sys- 
tematic 0.3 mag offset of the Clewley et al. My (BHB) 
relation with respect to the globular clusters. Therefore, 
we adjust the zero-point 0.3 mag brighter: 

Mv{BHB) = 1.252-(-0.214[Fe/H]-4.423(B-T/)o+17.74(B 

Although our zero-point adjustment suggests that the 
error in the absolute scale of BHB luminosities may be 
substantial, we expect that the relative BHB luminosities 
are precise to better than 10% for our sample. 

There are also 526 stars with early A spectral types be- 
tween B8 and A3 that are not BHB stars (see Figure [7]). 
These high surface-gravity stars likely comprise a mix of 
old blue stragglers and young main-sequence stars: two 
thirds of the early A-type stars have low mean metallicity 
[Fe/H] ~ —0.9; one third are consistent with s olar metal- 
licity [Fe/ H]=0. We use the My jA) relation of lSaraiedinil 
([l993§ and lKinman et al] (119941 ) to estimate luminosities 
for the 526 early A-type stars and the 124 BHB/A stars: 

Mv{A) = 1.32 + 4.05(S-y)o - 0.45[Fe/H] (2) 

This relation is based on a fit to globular cluster blue 
stragglers of similar spectral type. 

We estimate distances using the calculated luminosi- 
ties and the observed magnitudes. We convert 2MASS 
J magnitudes to V magnitudes by taking our [B — V)q 
estimate and lo oking up the corresponding iV — J)o in 
iKenvon fc Har tmanu ( 1995) for a star of that color. This 
conversion adds additional uncertainty to our distance 
estimates. Thus a typical BHB star with [Fe/H] = —1.7 
and {B-V)a = 0.07 has My (BHB) = +0.65 and a dis- 
tance of 9.3 kpc at our limiting magnitude J — 15.5, with 
typical uncertainty of 9% in distance. By comparison, 
a non-BHB, early A-type star with [Fe/H] — —0.9 and 
{B-V)a = 0.07 has My (A) = +2.0 and a distance of 
5.0 kpc at our limiting magnitude, with a typical uncer- 
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Fig. 9. — Mean rotation {upper panel) and line-of-sight velocity 
dispersion (lower panel) of 1225 BHB and early A stars located 
8 < i? < 11 kpc [solid lines). We divide the sample into thirds 
by [Fe/H] and find that the most metal-rich third [Fe/H] > -0.69 
{dotted lines) has systematically higher rotation velocity and lower 
velocity dispersion than the most metal-poor third [Fe/H] < —1.63 
{dashed lines). 



-]f)i^nt3-5cif3l(2f^^iSfJ)^.istance. These distance uncertainties 
do not include systematic errors. 

Stars of later spectral type than the BHB / early A 
stars are intrinsically less luminous objects in the nearby 
disk; we exclude these objects from our analysis and do 
not calculate their luminosities and distances. 

3.6. Mean Galactic Rotation 

Previous surveys provide conflicting measurements 
of the stellar halo r otation: it may be prograde 

fiiba fc BeersI \200d: 'Sirk o et all l2004bl). retrograde 
a i cwsk i ll99a~lM ai cwsk i et all 119961: ICarnev et all 
'1996; Soagna et alj i2003i : iKinman et aP '2004', '20071 
CaroU o et al. 2007) , or noth in g at all (Lavden et 
1996 :iGould fc Popowski 1998: Marti n fc Morrisonlll998f 
IGilmore et al.ll2002 : iBrown et al. . ,2005f ). Curiouslv. mea- 
surements of retrograde rotation mostly come from sur- 
veys of the north Galactic pole. Differences in the ob- 
served halo rotation may also arise from the manner in 
which different samples have s elected from inner-ha lo and 
outer-halo populations, which [CaroUo et al] ()2007f ) argue 
have quite different rotation characteristics. Our large 
area survey is ideal for measuring the bulk rotation of 
the inner stellar halo. Our stars are selected photometri- 
cally without any kinematic bias, and our survey is 100% 
complete over two large contiguous regions. 

We calculate mean , rotation using the formalism of 
iFrenk fc Whit3 (fl980h . We assume the stars are in pure 
rotation with uniform velocity about the rotation axis 
of the Galactic disk, and project the observed radial 
velocity onto t h e azim uthal (rotation) direction. The 
iFrenk fc White! (|1980f ) formalism provides an estimate 
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of both rotation velocity Vrot and the line-of-sight dis- 
persion aios for a set of stars. 

Figure [5] plots the results of this analysis as a func- 
tion of \Z\. We bin in \Z\ by binning together 78 stars 
ordered in \Z\, and moving through the sample in steps 
of 10 stars. We adopt this approach to avoid any ef- 
fects of arbitrary placement of bins on the results; bins 
are typically ^0.6 kpc in size, growing to > 1 kpc at 
|Z| > 5 kpc. We consider three different samples of stars: 
the clean sample of 655 BHB stars (solid line), a com- 
bined sample of 344 metal-poor BHB and early A-type 
stars with [Fe/H] < — 1.8 (dashed hue), and a combined 
sample of 541 intermediate metallicity BHB and early 
A-type stars with -1.8 < [Fe/H] < -0.6 (dotted line). 
The latter two metallicity cuts are intended to select halo 
and thick disk stellar populations, respectively, following 
iChibafc BeCTi(l2000l) . 

All stars rotate well below the solar value; Vrot drops 
monotonically with \Z\. This conclusion is valid for stars 
in the region 8 < i? < 12 kpc and 2 < |Z[ < 9 kpc. A 
linear least squares fit to the clean sample of BHB stars 
yields Vrot(BHB) = (-28 ± 3.4)|Z| -I- (175 ± 16) km s^^. 
The observed velocity gradient is statistically identical 
to the —30 ± 3 km s~^ gradient fo und for thick disk stars 
in the region < \Z \ < 4 kpc bvl Chiba fc Beers! (|2000D 
and by Girar d et all (j2006) (see Fig. (H). What is re- 
markable, however, is that all three of our samples show 
the same velocity gradient. The metal-poor stars have 
only a marginally shallower velocity gradient and lower 
zero-point Vrot{metal poor) = ("20 ± 3)|Z| -(- (129 ± 13) 
km s~^ than the intermediate-metallicity stars with 

Vrot{intermediate metallicity) — ( — 24±4) | Z | + ( 166± 13) km 

s~^. Thus the mean rotation velocities suggest that our 
samples contain significant numbers of thick disk stars 
with \Z\ < 5 kpc. 

The BHB and A stars located at 1.5 < \Z\ < 3 kpc 
have Vrot ~ 100 km s~^, consistent with the rotation 
measured from faint F and G stars at similar distances 
Gilmore e t al. 2002; Wyse et al. 2006). Gilmorc ct afl 



-0.5 



2002t l argue that a single, coherent thick disk should 



have a constant rotation velocity far from the plane, and 
thus the observed intermediate ~100 km s~^ rotation is 
evidence for a merger origin for the thick disk. How- 
ever, our sample clearly contains a mix of stellar popu- 
lations with different kinematics: the metal-poor sample 
has a systematically larger line-of-sight velocity disper- 
sion than the intermediate metallicity sample at a given 
\Z\. 

The lower panel of Figure [91 plots the line-of-sight 
velocity dispersion aios of the stars as a function of 
\Z\. The line-of-sight velocity dispersion of the clean 
BHB sample is statistically consistent with a constant 
value of 103 ± 6 km s""'^. This is consistent with the 
~100 km s'~^ velocity dispersion measured for thick disk 
stars at the sa me depth towards the south Galactic pole 
(jGirard et al.| [2006l. Yet the metal-poor sample of stars 
has aios = -1-117 ±10, consistent with a more halo- 
dominated population. Clearly, our BHB stars are a 
mix of thick-disk and halo populations. We note that 
the drop in velocity dispersion at \Z\ > 4 kpc seen in 
the metal-poor stars, though not statistically significant, 
may in fact be a real feature due to velocity structure in 
the sample. 




Fig. 10. — Position and metallicity of our clean sample of BHB 
stars. Solid squares are the average metallicity found in bins of 
50 stars. A linear least squares fit of the thick disk metallicity 
gradient 2 < \Z\ < 4 kpc finds d[Fe/H]/dZ = -0.03 ± 0.05. The 
mean metallicity of BHB stars with 5 < |Z| < 9 kpc is [Fe/H]= 2.0. 



We find no evidence for significant rotation of the inner 
halo in the range 5.5 < \Z\ < 9 kpc. Figure [9] shows 
that the formal Vrot values dip below zero in the region 
6 < |Z| < 8 kpc, but the uncertainties are large (note 
the error bars in Fig. [9]). Interestingly, the final bins 
of both the clean BHB sample and the [Fe/H] < —1.8 
BHB and A star sample have values around zero. In 
the volume 5.5 < |Z| < 9 kpc, the clean BHB sample 
contains 131 stars with a formal Vrot — —4 ± 31 km s~^, 
while the metal-poor sample contains 76 stars with a 
formal Vrot = —3 ± 37 km s"^. 

3.7. The Metal-Weak Thick Disk 

The thick disk is generally understood to have a metal- 
licity distribut ion th at peaks around [Fe/H] ^ — 0.7 (e.g. 
iGilmore et all 119951: lAllende Prieto et al.l |2006[ ). This 
conclusion is at odds with the low metallicity of our BHB 
stars located at |2'| < 4. Figure [TUl plots the position Z 
and metallicity of our clean sample of BHB stars. Solid 
squares are the average metallicity found in bins of 50 
stars. BHB stars located at 2 < |Z| < 4 kpc have an av- 
erage metallicity of [Fe/H] = —1.7; including the BHB/ A 
stars increases the average metallicity only slightly, to 
[Fe/H] = —1.4. The average metallicity of BHB stars 
in our sample with 5 < |Z| < 9 kpc is [Fe/H] = -2.0. 
The small numbers of stars in this region prevent a de- 
termination of a significant metallicity gradient, however 
it is interesting that the average metallicity is lower than 
one might have expected from a canonical halo popula- 
tion with peak metallicity [Fe/H] — —1.6 (jCarnev et al.l 
Il996f ). It may be that we are seeing evidence for the tran- 
sition of inner- t o outer-halo populat ions in this interval, 
as suggested bv lCarollo et all (|2007D . 

One explanation for both the low metallicity and the 
thick-disk- like kinematics (Figure [51) of our BHB stars at 
\Z\ < 4 kpc is p rovided by the existence of a metal-weak 
thick disk fe.g.. iNorris et al.lll985l: [Morrison et al.lll990l : 
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IChiba k Be"CTil200d : [Beers et al.ll2002[ ). BHB stars are 
associated with metal-poor populations, such as globular 
clusters, thus it may be that our survey for BHB stars 
preferentially traces the metal-weak thick disk. In any 
case, the existence of very metal-poor BHB stars with 
thick-disk kinematics presents another clue for formation 
scenarios of the Milky Way. 

A linear least-squares fit to the BHB stars located 
at 2 < |Z| < 4 kpc reveals a weak metallicity gradi- 
ent, d[Fe/H]/dZ = -0.03 ± 0.05 (sohd hne in Figure 
nop , consistent with zero. Previous studies find no ev- 
idence for a vertical metallicity gradien t in the more 
commonly studied (metal-rich) t hick disk l|Gilmore et al.l 
fl995t lAllende Prieto et aD[200l . 

4. LUMINOSITY FUNCTION OF BHB STARS 

Understanding the luminosity function of field BHB 
stars is important for interpreting our maps of the Galac- 
tic halo. BHB stars have a distribution of intrinsic lumi- 
nosities, thus we sample different luminosity BHB stars 
to different depths. The luminosity function describes 
the number of BHB stars per unit volume in the lumi- 
nosity interval My to My -|- dM. While we must in- 
fer a star's luminosity from its color and metallicity, the 
Mv{BHB) relation (Eqn. [T]) specifies only how a partic- 
ular color and metallicity map to a particular My. We 
emphasize that it is the observed distribution of BHB 
colors and metallicities that determines the form of the 
luminosity function. Paper H presents a more extensive 
discussion of this issue and the underlying physics. 

We calculate the fi e ld BH B luminosity function using 
the lEfstathiou et al.l (jlQSST ) non-parametric maximum- 
likelihood method. An important feature of this method 
is that the density terms drop out, thus the luminosity 
function calculation is unbiased by density variations. 
In other words, the maximum- likelihood method does 
not require knowledge of the halo density distribution 
p{R, Z), it only requires that the luminosity function is 
independent of position in the volume sampled. Because 
stellar density varies with position in the Milky Way, we 
compute only the form of the luminosity function and 
arbitrarily normalize the luminosity function to one. 

Figure [Tl] plots the luminosity function of the clean 
sample of BHB stars (dashed line). The luminosity 
function rises steeply at bright luminosities, peaks at 
A/y = 0.64, and falls rapidly with a tail at faint luminosi- 
ties. Although our statistics are greatly improved over 
the field BHB luminosity function measured in Paper II, 
we caution that our BHB sample is incomplete for stars 
{J-H)o>0.1. 

In principle, we can correct for our sample incom- 
pleteness. The Paper I sample is complete over a much 
broader range of color than our 2MASS-selected sample. 
Thus the Paper I sample can provide us with the dis- 
tribution of BHB colors with (J — H)o > 0.1 that are 
missing from our 2MASS-selected sample. We estimate 
the luminosities of the missing, redder BHB stars as fol- 
lows. First, we determine the distribution of BVO colors 
of our stars as a function of (J — H)o. Second, we deter- 
mine the distribution of [Fe/H] for BHB stars with colors 
near {J — H)o = 0.1. Third, we construct cumulative dis- 
tributions of BVO and [Fc/H] from our observations, and 
then sample these distributions to obtain the expected 
distribution of My's for the missing stars. Finally, we 



(X 



CO 



' ' ' I ' ' 

BHB observed 
BHB corrected 
including BHB/A 




r-1 1 i~i 



0.5 1 
My [mag] 



1.5 



Fig. 11. — Luminosity function of field BHB stars, calculated for 
the observed sample (dashed line) and corrected for completeness 
(solid line). If we include the BHB/A stars with [Fe/H] > -0.6 
(dotted line) a faint tail appears. However, it is clear that hot, 
intrinsically faint extended BHB stars are not a significant fraction 
of the field BHB population. 



correct the luminosity function bins for the appropriate 
fraction of missing stars as determined from the Paper I 
sample. 

Figure [TT] plots the BHB luminosity function corrected 
for incompleteness (solid line) . The effect of the incom- 
pleteness correction is to increase the fraction of redder, 
more luminous stars, and thus shift the peak of the dis- 
tribution to My — 0.60. We also plot the corrected lu- 
minosity function for the combined sample of BHB and 
BHB/A stars. Because the BHB/A stars have [Fc/H] 
> —0.6 and bluer colors, on average, than the BHB stars, 
they are intrinsically under- luminous and fall entirely in 
the faint tail of the luminosity function. These stars are 
possibly hot extended BHB stars, though such a strong 
preference for high metallicities is not observed in globu- 
lar clus ters. For example, NGC 7078 has [Fc/H] = -2.25 
(jHarriSi 1996) and contains a large number of extended 
BHB stars. Even if the BHB/A stars are all BHB stars, 
it is clear from Figure [TT] that extended BHB stars are 
not a significant fraction of the field BHB population. 

5. DENSITY DISTRIBUTION OF FIELD BHB STARS 

Our sample of BHB stars can potentially provide an ex- 
cellent measure of the density distribution of the thick- 
disk and inner-halo stellar populations. Traditionally, 
the Galaxy's density distribution is measured with star 
counts. The star-count technique is powerful because can 
use photometric catalogs containing millions o f stars (e.g. 
ISiegel et al.ll2002l : iLarsen k HumpIirev£ll2003D . However, 
stellar populations are a complex function of both color 
and apparent magnitude. Thus star-count techniques 
suffer from uncertainties in stellar color-luminosity re- 
lations, as well as contamination from binaries and non- 
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Fig. 12. — a) Best-fit density distribution (solid line) and the observed density of BHB stars at R = 9 kpc (points). The extrapolation 
(dashed line) suggests there are 78 ± 30 kpc"'' halo and thick disk BHB stars at (R, Z) = (9, 0) kpc. b) The fraction of thick disk and halo 
13HB stars in our survey volume located 8 < iJ < 10 kpc. 



stellar objects. By comparison, our survey provides a 
very clean sample of spectroscopically identified BHB 
stars. Although the numbers of BHB stars is much 
smaller than samples of stars used by star counts, our 
spectra provide precise metallicity and distance determi- 
nations for every star. 

We begin by considering the volume of space sampled 
by our survey. Using the BHB luminosity function, we 
calculate the fraction of stars at a given distance that 
fall within our survey magnitude limits. We expect that 
our survey is more than 50% complete for BHB stars in 
the range 2.5 < d < 9.5 kpc. This range of heliocen- 
tric distance samples the region 8 < i? < 12 kpc and 
2 < \Z\ < 8.5 kpc for our predominantly high Galactic- 
latitude survey region. Thus we restrict our density dis- 
tribution analysis to the above ranges of d, R, and Z. 
There are 544 BHB stars that fall within these ranges. 

Before calculating the density distribution, we correct 
the observed BHB sample for incompleteness. First, the 
BHB luminosity function tells us the fraction of stars 
missing at each distance d. Because of our restriction in 
distance, this correction applies to only a handful of stars 
near the boundaries of the sample. Second, the ratio of 
the observed and corrected luminosity functions tells us 
the fraction of BHB stars missing at each My because 
of our color selection; we weight stars at a given My 
appropriately. Finally, we use Figure [3] to estimate the 
fraction of faint BHB stars missing because of increased 
photometric errors, and give additional weight to stars 
with J > 14.5. 

We assume the density distribution is a sum of a thick- 
disk and halo population with the canonical forms: 

p{R, Z) = po,thick exp (-Z/hz) exp {-R/hfi)+po,haio/ (oo 

(3) 

where Rg = \J R? ^ [z c/aY and c/a is th e halo axial 
ratio. These are the same relations used by ISiegel et al.l 
()2002| ) in their fits to star counts. Selecting a power- 
law halo instead of a de Vaucouleurs halo is mostly a 
cosmetic choice; our BHB sample provides very little 
leverage on the halo density profile. Our sample is a 
high Galactic-latitude sample best suited to measuring 



hz and the relative normalizations of thick disk and halo 
B HB star dens i ties. 

iGirard et al.l (|2006D caution that distance uncertain- 
ties, convolved with the sharply falling density distribu- 
tion of stars, can alter the "observed" density distribu- 
tion from the actual, intrinsic form. Thus, we mimic their 
procedure, and artificially partition each star into 100 
positions in distance, with a distribution of distances de- 
scribed by a Gaussian distribution around the best value 
for each star. We then apply our limits in d, R, and Z to 
the subunits. This procedure allows stars that would oth- 
erwise be excluded by distance cuts to contribute a small 
amount of weight appropriate to the uncertainty in their 
distance estimate. We bin the subunits into volumes at 
fixed intervals R and Z, and perform our de nsity fits to 
these bins using minimization techniques ([Press et al.l 

mn). 

We start by testing fits to the different components 
of our BHB sample. The density distribution of stars 
located at —4 < Z < —2 kpc and 2 < Z < 4 kpc are 
very similar. In the final fit we consider stars above and 
below the plane together as a function of \Z\. We try 
different bin sizes and different ranges of R and Z, and 
find that stars located at 2 < |Z| < 4 kpc prefer thick- 
disk scale lengths in the ranges 2.5 < /ij^ < 4 kpc and 
\ < hz < 2. Unfortunately, stars located farther out, 
at 5 < \Z\ < 8.5 kpc, provide very little constraint on 
the form of the halo profile. If we fix the core radius 
to oo = 6.3 kpc (Girard et all 20061) and the halo axia l 
ratio to c/a = 0.7 (iRobin et aLll2000! : 'Siegel et al."200^, 
our sample prefers an halo power-law index in the range 
2.5 < n < 3. 

We fit Equation [3] to the full BHB sample, holding 
'^gc}/a ~ 0.7 fixed and fitting the other 6 parameters. Table 
[T] gives the best-fit parameters and Figure [T^ shows the 
results. Figure [T2k plots the observed density of BHB 
stars at i? = 9 kpc, the best-fit density distribution (solid 
line), and the thick-disk and halo components (dotted 
lines) . Figure [l2b plots the fraction of thick-disk and 
halo stars in our survey volume located 8 < i? < 10; the 
components contribute equal fractions of BHB stars at 
\Z\ - 4.5 kpc. 
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TABLE 1 
BHB DENSITY 
DISTRIBUTION 



Par am 


Value 


Units 


PO,thick 


960 ± 170 


kpc~^ 


hz 


1.26 ±0.1 


kpc 


ha 


3.5 ±0.5 


kpc 


PO,halo 


1040 ± 180 


kpc~^ 


ao 


5.8 ±3 


kpc 


n 


-2.5 ±0.5 




c/a 


= 0.7 





We map out contours of to understand the uncer- 
tainty in our best-fit parameters. We caution that our 
sample provides little constraint on oq, and that there 
is a significant degeneracy between the thick-disk scale 
length, hn, and the halo power law index n. Star-count 
models give si milar results to our bes t-fit halo pow er-law 
index n ~ 2 5 (iRobin et al.| r2000: Sic gel et al.ll2002l e.g.). 
In contrast. lChiba fc Beers! (j2000.') find n ~ 3.5 from their 
spectroscopic sample of metal-poor stars. Given the un- 
certainties, however, our BHB stars provide no significant 
constraint on the shape of the halo. What we can mea- 
sure with certainty is the normalization of halo to thick 
disk BHB stars. 

The density of halo and thick disk BHB stars is 
104±37kpc-3 near the Sun (R,Z) = (8,0 ) kpc, in good 
agreement with iGreen fc Morrison (I1993D 's lower limit 
of 51 ± 17 kpc""^. iKinman et al.l ( 19941 ) report a three 
times smaller density of 30 kpc~'^, but their sample has 
little constraint on thick disk BHB stars that dominate 
the BHB density near the Sun. We find that the relative 
normalization of halo to thick disk BHB stars is 4 ± 1% 
near the Sun. 

Our high-latitude sample also provides a good con- 
straint on the vertical density distribution of BHB stars 
in the thick disk. To obtain a self-consistent picture of 
the relationship between the thick-disk scale lengths hz 
and hn, we fix the form of the halo power law (see Table 
[1]), vary hz and hfj across a grid of values, and fit only 
the normalizations. This approach results in contours of 
illustrated in Figure 1131 The contours do not corre- 
spond to exact significance levels, but we have chosen the 
inner contour to match our best estimate of la signifi- 
cance based on boot-strap resampling. The asterisk in 
Figure [T^ marks our best-fit values of hz and hR. 

Previous star-count models for the thick disk find ei- 
ther a larg e scale height (1.2 - 1 . 4 k pc) and low nor- 
malization (iGilmore fc ReidI 1 1983t [Morrison et al.l l2000l: 
iReid fc Maiewskil 119931 : lJuric et al.l l2005f) or a~smaller 
scale height (0.75 - 1.0 kpc) and high normalization 
jRobin et al."1996'; 'Siegel et al."2002'; 'Robin et al .ll200l 
Cabrcra-Lavcrs et al. 2005; Du et al. 2006; Girar d et al.l 
20061) . Our scale height hz = 1.26 ± 0.1 is consistent 



with the larger scale heights. 

The thick-disk scale length is interesting because 
it determines the rotational-equilibrium of thick disk 
stars. Star-coun t estim ates r ange from 2.5 ± 0.3 
kpc (i Robin et al.' 1996, ' 2003f) to 4.3 ± 07 kpc 
(jLarsen fc Humphreva .200311 . Our scale length, hr = 
3.5 ± 0.5 kpc, falls in the middle of this range, similar to 




h„ [kpc] 

Fig. 13. — contours for thick disk scale lengths; the asterisk 
marks our best-fit values for hz and hji. We calculate these con- 
tours by fixing the the halo power law and fitting the thick disk and 
halo normalizations. The inner contour matches our best estimate 
of 1(T significance based on boot-strap resampling. 



determ inations from lSiegel et al.l (|2002f) and lJuric et al.l 
(|2005l ). This agreement shows the power of a clean, even 
if small, spectroscopic sample. If our BHB stars prefer- 
entially trace the metal- weak thick disk, then the density 
parameters in Table 1 reflect the density distribution of 
the metal-weak thick disk. 

6. CONCLUSIONS 

We discuss a complete spectroscopic survey of 2414 
2MASS-selected BHB candidates over 4300 deg^ of sky. 
We identify 655 BHB stars in this non-kinematically se- 
lected sample. The luminosity function of the field BHB 
stars has a median value of My = 0.65 and a small tail 
extending to My > 1, but shows very few extended BHB 
stars. 

The BHB stars located at |Z| < 4 kpc have a mean 
Galactic rotation and density distribution remarkably 
consistent with a metal-weak thick-disk population. The 
\Z\ < 4 kpc BHB stars have a mean metallicity of [Fc/H] 
= —1.7, a velocity gradient of dvrot/d\Z\ — —28 ± 3.4 
km s~^, and a vertical scale height of hz = 1.26 ±0.1 
kpc. We infer a space density of 100±36 kpc""^ thick disk 
BHB stars near the Sun. RR Lyrae stars, by comparison, 
have a muc h less prominent disk component near the Sun 
llMartin fc" Morri son 1998). The existence of metal-poor 
BHB stars with thick-disk kinematics and scale heights 
present another clue for formation scenarios of the Milky 
Way. 

The BHB stars located at 5 < \Z\ < 9 kpc have a 
mean Galactic rotation and density distribution consis- 
tent with a predominantly halo population. The halo 
BHB stars have a mean metallicity of [Fe/H] = —2.0, 
a mean Galactic rotation of —4 ± 31 km s~^, and rela- 
tive normalization of 4 ± 1% with respect to the thick 
disk near the Sun {R, Z) ~ (8, 0) kpc. This is the best 
currently available measurement of the relative normal- 
ization of the halo and thick disk, an important quantity 
that enters into all models of the local structure of the 
Galaxy. 

In the future, having established the global properties 
of our survey, we hope to analyze the BHB sample for 
structure in space, velocity, and metallicity. 
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TABLE A2 
PHOTOMETRY 
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Note. — 


Table EH is 


presented in its entirety in the electronic edition of the Astronomical Journal. A portion 
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TABLE A3 

SPECTROSCOPIC AND STELLAR PARAMETERS 
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the electronic edition of the Astronomical 



Journal. A portion is shown here for guidance and content. 
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APPENDIX 
DATA TABLES 

Tables IA2I and IA3I contain the photometric and spectroscopic measurements for the 2414 2MASS-selected BHB 
candidates. Our survey boundaries inc lude 100 objects previously observed as p art of P apers I and II. We include the 
previously published objects in Tables [X2l and [X3l for completeness. Tables E2] and E3l are presented in their entirety 
in the electronic edition of the Astronomical Journal. A portion of the tables are shown here for guidance regarding 
their format and content. 

Table IA2I summarizes photometry and positions. Column (1) is our identifier. The designation CHSS stands for 
Century Halo Star Survey and is chosen to be unique from previous surveys. Column (2) is the J2000 right ascension 
in hours, minutes, and seconds. Column (3) is the J2000 declination in degrees, arcminutes, and arcseconds. Column 
(4) is the extinction-corrected 2MASS Jo magnitude. Columns (5) and (6) are the ex tinction -corrected 2MASS colors 
(J - H)o and {H - K)q. Column (7) is the E{B - V) reddening value from Schlcge reTall (|T998). Columns (8) and 
(9) are the Galactic coordinates, in degrees. Column (10) is our BVO color predicted from 2MASS photometry and 
Balmer line strengths (and SDSS photometry, where available). 

Table R3l summarizes the spectroscopic and stellar parameters. We include all the DA white dwarfs and subdwarfs 
in this table, but we omit their stellar parameters as our analysis is meaningless for these objects. We also omit stellar 
parameters for a few dozen objects with unusually low signal-to-noise spectra. Column (1) is our identifier. Column 
(2) is the heliocentric radial velocity in km s^^. Column (3) is the spectral type, where B0=10, A0=20, F0=30, and 
so forth. Column (4) is the effective temperature in K. Column (5) is the surface gravity in cm s~^. Column (6) is 
the metallicity given as the logarithmic [Fe/H] ratio relative to the Sun. Column (7) is our classification: BHB = blue 
horizontal branch star, BHB/A = possible blue horizontal branch star with [Fe/H] > —0.6, A = high surface gravity, 
early A-type star, DA — DA white dwarf, sd = subdwarf. Column (8) is the absolute My magnitude estimated from 
Equations [Hand [2] for BHB and A stars, respectively. Column (9) is the estimated distance in kpc. Absolute magnitude 
and distance estimates are only provided for BHB and A-type stars, as described in Section 3.5. 



